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• of  generalized  coordinates  to  twelve.  Two  check  cases  are  presented  to  establish  the  accuracy 
of  the  finite  element. 


PREFACE 


This  work  is  part  of  an  ongoing  investigation  of  the  b< ' avior  of  pressurized  structural 
elements  under  load  with  the  ultimate  objective  of  making  x>ssib!e  the  use  of  such  elements 
in  the  support  structure  of  Army  tents.  The  work  WdS  funded  under  the  in-House 
Laboratory  Independent  Research  program  as  a work  unit  entitled,  "Study  of  the  Stability 
of  Pressure  Stabilized  Arches  and  their  Structural  Assemblies."'  In  the  reference  citations 
the  organizations  "US  Army  Natick  Laboratories"  and  "US  Army  Natick  Development 
Center"  refer  to  the  organization  now  called  the  "US  Arm/  Nnick  Research  and 
Development  Command." 
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PRESSURE  STABILIZED  BEAM  FINITE  ELEMENT 


INTRODUCTION 

A systems  analysis  of  Army  needs  for  shelters  in  the  1985  time  frame  revealed  that 
the  requirements  for  lightweight  tentage  of  low  packaged  bulk  with  minimum  setup  and 
disassembly  time  might  be  most  effectively  achieved  by  a tent  with  a frame  of  highly 
pressurized  (compared  with  present  air-supported  tents)  structural  elements  supporting  a 
lightweight  fabric  bamier,  as  illustrated  in  Figure  1.  Since  the  structural  elements  are 
pressure  stabilized,  they  can  be  fabricated  from  lightweight,  flexible  materials  such  as  cloth, 
thus  reducing  the  weight  and  bulk  of  the  transport  configuration.  The  use  of  bladders 
or  coated  fabric  for  these  elements  should  provide  sufficient  air  retention  capability  to 
eliminate  the  need  for  the  dedicated  air  supply  used  with  current  air-supported  tents. 

To  successfully  develop  this  concept,  a structural  design  capability  was  needed.  To 
meet  this  need  an  investigation  of  the  structural  behavior  of  pressure-stabilized  beams  and 
arches  was  carried  out  and  reported  in  references  1,  2,  and  3.  These  references  3ive 
the  development  of  theories  with  experimental  confirmation  to  predict  the  structural 
behavior  of  pressure  stabilized  beams  and  arches  as  a function  of  geometry,  pressure  level 
and  material  properties.  The  theories  give  adequate  predictions  of  both  the  deformation 
and  the  load-carrying  capability.  These  theories  thus  serve  to  characterize  the  behavior 
of  the  structural  elements  and  can  be  used  as  a guide  to  design. 

Since  these  structural  elements  are  to  be  used  in  structural  frames  for  tents,  a more 
complete  design  capability  is  needed;  that  is,  the  capability  to  predict  the  deformation 
and  load-carrying  capability  of  frame  assemblies  made  of  pressure-stabilized  elements  and 
loaded  through  a fabric  barrier.  Thus,  a structural  analysis  method  capable  of  treating 
the  interaction  of  assemblages  of  frame  elements  and  the  interaction  of  this  frame  with 
the  fabric  barrier  is  needed.  The  finite  element  procedure  is  well  suited  to  such  calculations, 
and  a computer  program  using  that  procedure  for  the  analysis  of  conventional  frame  tents 
was  available.  This  program,  which  is  described  in  references  4 and  5,  includes  a 


1.  Earl  C.  Steeves;  A Linear  Analysis  of  the  Deformation  of  Pressure  Stabilized  Beams;  US 
Army  Natick  Laboratories;  Technical  Repor  /5-47-AMEL,  1975  (AD  A006493). 

2.  Earl  C.  Steeves;  Behavior  of  Pressure  Stabilizeo  Beams  Under  Load;  US  Army  Natick 
Development  Center;  Technical  Report  75-82-AMEL;  1975  (AD  A010702). 

3.  Earl  C,  Steeves;  The  Structural  Behavior  of  Fressure  Stabilized  Arches;  US  Army 
Research  & Development  Commend;  Natick/TR— 78/018;  1978, 

4.  Paul  J,  Remington,  John  C,  O'Callahan  and  Richard  Madden;  Analysis  of  Stresses 
and  Deflections  in  Frame  Supported  Tents;  US  Aony  Natick  Laboratories;  Technical 
Report  75-31;;  1974  (AD  A002072). 

5.  Paul  J Remington,  John  C.  O'Callahan  and  Richard  Madden;' Finite  Element  Analysis 
of  Scale-Model  Frame-Supported  Tents;  US  Army  Natick  Research  & Development 
Command;  Technical  Report  76— 21— AMEL,  1975. 
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conventional  beam  element  and  a nonlinear  membrane  element  for  modeling  the  fabric 
barrier.  Since  a pressure  stabilized  beam  element  was  not  included,  it  was  decided  to 
develop  such  an  element  and  adapt  it  to  this  existing  computer  program.  The  purpose 
here  is  to  report  this  work.  This  will  include  the  development  of  the  finite  element, 
its  adaptation  to  the  computer  code  and  some  check  cases. 

DEVELOPMENT  OF  FINITE  ELEMENT 

We  proceed  here  with  the  derivation  of  the  finite  element  with  the  intended  result 
a 12  X 12  stiffness  matrix  associated  with  generalized  coordinate  set  comprised  of  the 
three  deflections  and  three  rotations  of  each  of  the  beam  end  points.  The  basis  for  this 
derivation  is  the  energy  principle  derived  in  reference  1.  This  energy  principle  which 
includes  bending  in  one  plane  and  axial  deformation  is  augmented  to  include  bending 
in  two  planes,  axial  deformation,  and  torsion.  The  displacement  fields  are  represented 
by  polynomials  in  the  coordinate  specifying  position  along  the  beam. 

Because  of  the  number  and  order  of  polynomials  used,  the  number  of  generalized 
coordinates  exceeds  12,  so  a condensation  is  carried  out. 

In  carrying  out  this  derivation  we  use  the  coordinate  system  and  sign  convention 
given  in  Figure  2.  Shown  are  the  sign  conventions  associated  with  the  basic  theory  and 
that  associated  with  the  finite  element.  The  potential  energy  expression  to  be  used  is 
that  given  in  reference  1 and  augmented  as  indicated  above  and  is  expressed  in  dimensional 
form  as: 


+ D{^): 

dx 


+ 


_ ,dv 
+ C< 

dx 


- UV»  - 

v)*  + C(_  + 0)* 

dx 


+ A ♦ p&='>  dx 

dx  dx  j 


This  expression  is  based  on  the  assumptions  of  uniaxial  strain  and  linear  expression  for 
the  shear  strain  in  terms  of  the  displacements.  Since  we  seek  a finite  element  with  which 
forces  will  be  applied  at  the  nodes,  the  distributed  forces  are  not  included  in  (1).  The 
potential  energy  expression  is  a quadratic  in  the  displacements  and  can  be  written  in  matrix 
form  as 
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FIGURE  2.  COORDINATE  SYSTEM  AND  SIGN  CONVENTION 


du  dup  dW_  ^ dd  dv  ^ ^ 

dx  dx  dx  dx  dx  dx 


C+P  c 


C C 


W d 0 > -j*  *2) 

1 1 dx  I 


C U 


where 


E = 27raC,  , 
C = ?raCi  3 
P = 7ra2P/2 


D = ^a^Ct  i 


G = 7ra3Ci 


In  these  definitions  C,  ( , C3J,  and  p are  respectively  the  elastic  modulus,  the  shear  modulus, 
and  the  pressure.  It  can  easily  be  seen  from  this  espression  that  the  axial,  torsional, 
and  two  bending  energies  are  uncoupled  and  thus  can  be  treated  independently 
informulating  the  finite  element  stiffness  matrix. 

du  , 

The  energy  resulting  from  axial  deformation  is  given  by  the  expression  which 

is  the  classical  expression.  With  the  use  of  a linear  representation  for  the  axial  displacement 
field  in  terms  of  the  end  point  axial  displacements,  U,  and  U2,  the  following  energy 
expression  is  obtained  upon  integration: 


2]  ^ 


E/L  — E/L 1 U, 


72  [_— E/L  E/L  j [U2J 

This  is  the  classical  stiffness  matrix  for  linear  axial  behavior. 


Similarly,  the  torsional  energy  is  given  by  G(£^)2  and  a representation  of  0 linear 
in  the  x coo, di.. ate  in  terms  of  the  end  point  rotation  0X  and  02(  yields  for  the  energy 


9 


G/L  V t 


0, 


G/L 

-G/L 


G/L 


(K 


(4) 


Examination  of  (2)  reveals  tnat  the  two  bending  problems  differ  only  in  the  algebraic 
signs  of  the  off  diagonal  elements.  This  being  so,  we  can  derive  the  stiffness  matrix 
associated  with  one  of  the  bending  problems  and  dednce  the  other  form  it  We  will 
treat  the  problem  described  by  W and  o The  formulation  of  the  bending  stiffness  matrix 
is  not  as  direct  or  simple  as  the  above  problems  or  the  classical  bending  problem  This 
results  from  the  use  of  two  independent  parameters,  W and  Q,  to  descriDe  the  bending 
problem.  The  use  of  four  generalized  coordinates  is  generally  sufficient  to  approximately 
describe  the  beam  bending  problem  and  is  though  to  be  adequate  here  also.  However, 
dividing  these  four  coordinates  between  the  two  independent  variables,  W and  q,  only 
allows  for  the  use  of  linear  aporcximaoon  of  these  variables'  w'thin  the  context  of 
polynomial  approximation.  That  is  not  believed  to  be  an  adequate  approximation,,  since 
for  W it  only  allows  description  of  rigid  body  displacements.  It  was,  therefore,  decided 
to  use  two  third  order  polynomials,  one  to  represent.  W and  the  other  to  represent  q 
This  gives  a stiffness  matrix  in  terms  of  eight  generalized  coordinates,  so  a condensation 
is  performed  to  reduce  the  numu_'  io  four.  The  bending  energy  is 


and  the  polynomial  approximations  are 
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In  equation  (6) 


(7) 


h,  = 1 - 3x2/L2  + 2x3/L3  h'j  = 6x2/L’  - 6x/L2 

h;  = 3x2/L2  - 2xVL’  h = 1 - 4x/L  + 3x:/L: 

h3  = x — 2x2/L  + X‘VL:  = 3x2/L2  — 2x/L 

h4  = xVL2  — x2/L 


The  elements  of  the  vector  of  order  eight  in  equation  (6)  are  the  transverse  deflection 
and  cross-section  rotation  and  their  derivations  evaluated  at  x = 0 and  x = L.  The 
subscript  1 denotes  evaluation  a*  x = 0,  and  suoscript  2 denotes  evaluation  at  x = L.: 
The  prime  denotes  the  derivative  of  the  variable.  Substitution  of  (6)  into  (5)  and 
integrafion  yields  the  following  energy  expression: 


Fa! 
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the  matrices  K,  , 

K 1 2 

(9) 


an  eighth-order  stiffness  matrix  which  we  would  like  to  reduce  to  fourth  order  to  make 
it  compatible  with  other  finite  elements  for  beam  bending  but  particularly  for  compatibility 
with  the  finite  element  described  in  references  4 and  5.  To  accomplish  this  the  set  of 
generalized  coordinates  contained  in  A,  are  chosen  to  be  the  set  to  be  used  and  those 
in  A2  the  set  to  be  removed.  It  is  also  assumed  that  the  applied  forces  do  work  only 
through  the  generalized  coordinates  A,  so  that  the  potential  energy  including  the  external 
forces  can  be  written 
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In  this  situation  the  following  transformation  of  coordina^s  can  be  used  to  carry  out 
the  desired  reduction  or  condensation 

V " I ](a,' 

= j 

A,  | -K2\  K,  2j  (11) 

with  I denoting  the  identity  matrix.  Substitjti  'n  of  (11)  into  (10)  gives  for  the  energy 
expression 

y2  A^SA,  - a\z  (12) 

Where  S is  the  desired  4x4  stiffness  matrix  expressed  in  terms  of  the  original  suomatrices 
as 

S = K,  , - Ki2  K;1,  k!2  (13) 

This  completes  the  development  of  the  stiffness  matrices  for  axial,  tort'onal,  and  bending 
deformation  and  the  complete  12  x 12  stiffness  matrix  associated  with  the  end  point 
displacements  and  rotations  can  be  assembled.  In  doing  so,  attention  is  paid  to  the  order 
of  the  generalized  coordinates  or  end  point  displacements  used  in  the  NONFESA  (Nonlinear 
Finite  Element  Structural  Analysis)  program  in  which  this  finite  element  will  be  used. 
The  desired  stiffness  matrix  appears  as  follows  in  the  expression  for  the  potential  energy: 
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I Ilf  i'll  1 1 k 'I its  s()  111  (14)  .iri>  Itic  rlumiints  ol  the  matrix  S given  in  M3)  In  assembling 

lilts  iinitnx,  count  I ws  been  taken  ol  the  algebraic  sign  difference  of  the  off  diagonal 
elements  in  the  bending  problem  associated  with  V and  i p that  was  discussed  previously 
in  eoniwtion  with  expression  (?)  The  matrix  in  (14)  is  the  desired  stiffness  matrix, 
,iihI  Ii.iviihi  il  we  ran  pioreed  with  its  adaptation  to  the  NONFfcSA 

In  (aityitig  out  this  adaptation,  two  options  were  open,  the.  first  being  to  create 
an  enlneiy  new  element  to  be  included  along  with  the  beam,  membrane  and  truss  elements 
alieady  in  the  piogiam,  and  the  second  being  to  modify  the  present  beam  element  by 
making  the  piessure  stabilized  beam  clement  an  option  on  the  conventional  beam  element. 
Because  ol  the  great  similarity  of  this  beam  element  to  the  conventional  one,  this  second 
option  was  chosen.  To  accomplish  this  adaptation  changes  were  required  in  both  the 
piogtom  and  the  input.  Program  changes  were  required  in  the  computation  of  the  stiffness 
matrix  and  in  the  computation  of  the  stresses  for  output.  The  beam  stiffness  matrix 
is  computed  in  the  subroutine  named  NEWBM,  and  in  this  subroutine  an  option  was  added 
for  calculation  of  the  stiffness  matrix  according  to  (14).  This  option  is  controlled  by 
the  parameter  JPRES.  For  the  calculation  of  the  stiffness  matrix  according  to  (14)  two 
additional  subroutines  are  needed,  ELB  to  calculate  the  stiffness  matrix  S and  GJSOL 
to  solve  a system  of  equations  by  Gauss-Jordan  reduction.  In  ELB  the  order  eight  stiffness 
matrix  of  (8)  is  calculated  according  to  the  formulas  for  the  elements  given  in  the  appendix. 
The  coi  iensation  of  the  matrix  is  carried  out  using  subroutine  GJSOL,  and  the  order 
four  matrix  is  returned  to  subroutine  NEWBM  for  assembly  into  the  12  x 12  stiffness 
matrix  along  with  the  other  elements  which  are  calculated  directly  into  the  matrix.  The 
reader  may  question  the  need  for  inclusion  of  the  subroutine  GJSOL  for  Gauss-Jordan 
reduction  when  one  is  already  present  in  NON  F ESA  for  solution  of  the  large  system  of 
equations  generated  by  the  program.  The  subroutine  already  included  is  especially  designed 
for  large  systems  of  equations  stored  on  tape  in  block  format  so  its  use  on  a fourth 
order  system  seemed  inefficient  and  overly  complex.  Thus  a very  simple  25  line  subroutine 
was  added  to  carry  out  the  condensation.  The  other  programming  changes  are  in  subroutine 
BEAM  where  the  stresses  are  computed  for  output  from  the  nodal  forces.  Because  the 
pressure-stabilized  beam  is  a circular  cylindrical  membrane,  membrane-type  stresses  are 
computed  as  follows: 

The  uniform  stress  resulting  from  axial  deformation  is  computed  by  dividing 
the  axial  nodal  forces  by  the  circumference  of  the  cross-section. 

The  maximum  and  minimum  stresses  rulting  from  bending  deformation  are 
computed  for  each  of  the  two  planes  of  bending.  The  magnitudes  of  these 
parameters  are  the  quantities  and  their  negatives  obtained  by  dividing  the  nodal 
moment  for  each  of  the  planes  of  bending  by  the  area  inclosed  by  the 
cross-section. 

The  average  stresses  from  shear  deformation  are  computed  by  dividing  nodal 
shear  force  associated  with  each  of  the  two  planes  of  bending  by  the 
circumference  of  the  cross-section.. 
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These  stiess  calculations  are  carried  out  for  both  ends  of  the  beam. 

The  only  input  that  changes  is  that  associated  with  the  beam  element,  with  these 
changes  being  on  the  material  property  card,  the  element  property  card  and  the  beam 
data  cards  of  reference  6.  For  the  pressure  stab,  zed  beam  the  needed  input  in  addition 
to  the  control  parameter  is  the  beam  length,  the  clastic  modulus,  Poisson's  Ratio,  the 
cross  section  radius,  and  the  inflation  pressure.  The  beam  length  is  computed  from  the 
beam  end  point  nodal  coordinates  and  this  does  not  constitute  a change.  The  elastic 
modulus  and  Poisson's  Ratio  are  entered  with  the  material  property  card,  with  the  elastic 
modulus  being  the  first  entry  following  the  identification  number  and  Poisson's  Ratio 
the  second.  The  elastic  modulus  required  is  of  the  membrane  type  having  the  units  of 
force  per  unit  length.  If  Poisson's  Ration  is  known,  it  can  be  entered,  but  typically 
the  parameter  that  is  known  is  the  membrane  shear  modulus.  In  the  program  the  shear 
modulus  is  computed  from  the  elastic  modulus  and  Poisson's  Ratio  according  to  the  relation 

C,3  = C, ,/2(1  + n)  (15) 

Thus,  when  the  shear  modulus  and  the  elastic  modulus  are  known,  the  Poisson's  Ratio 
to  be  entered  as  the  second  entry  on  the  material  property  card  should  be  calculated 
by  the  expression  by  solving  (15)  for  n in  terms  of  C,  , and  C,3.  The  element  propeiy 
card  contains  an  identification  number  followed  by  eight  numerical  parameters.  For  the 
pressure-stabilized  beam  elements,  the  first  two  of  the  numerical  parameters  are  the 
cross-section  radius  and  the  inflation  pressure,  in  that  order.  The  fifth  of  these  parameters 
should  be  entered  as  a negative  number  as  this  is  used  to  control  the  option  on  the 
computation  of  the  stresses  from  the  nodal  forces.  The  remaining  five  parameters  are 
dummies  for  the  pressure-stabilized  b'am  element,  and  nonvamshing  entries  should  be  made 
for  them.  Control  of  the  pressure-stabilized  beam  option  is  accomplished  by  the  parameter 
JPRES  on  the  beam  data  card  (there  is  one  of  these  cards  for  every  beam)  If  a given 
beam  is  to  use  the  pressure-stabilized  beam  option  the  PJRES  must  be  entered  as  unity 
in  column  68  of  the  beam  data  card,  otherwise  it  is  to  be  left  blank  or  entered  as  zero. 
If  the  pressure-stabilized  beam  element  is  exercised  then  the  material  property  and  element 
property  numbers  entered  on  the  beam  data  card  must  correspond  to  the  identification 
numbers  corresponding  to  material  and  element  property  sets  for  pressure-stabilized  beams 
as  described  above. 


6 John  C.  O'Callahan  NONFESA  - Nonlinear  Finite  Element  Structural  Analysis  Code 
for  the  Analysis  of  Stresses  and  Deflections  in  Frame  Supported  Tents,  Bolt  Beranek 
and  Newman,,  Inc.,  Report  No.  2803,  1975. 
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CHECK  CASES 


To  establish  ilio  validity  of  the  pressure  stabilized  beam  finite  element  two  check 
cast's  are  piesented  line  A straight  beam  loaded  a ns  midpoint  in  two  perpendicular 
planes  and  a eurulai  aicli  loaded  in  its  plane  are  the  check  cases.  The  finite  element 
calculations  aie  compared  with  the  exact  solutions  of  the  theory  given  in  references 
4 and  b 

Tne  simplest  structure  th.it  can  be  modeled  using  the  pressure  stabilized  team  finite 
element  is  a straight  beam,  as  shown  in  Figure  3.  We  have  taken  a beam  i20  cm  long 
having  a circular  cross  section  of  3 cm  radius  and  divided  it  into  six  elements  of  equal 
length.  The  element  numbers  h Figure  3 are  inclosed  in  the  square  symbols  and  the 
nodes  are  numbered  one  through  8 and  denoted  by  the  solid  circular  symbols.  Node 
eight  is  a dummy  node  require  to  establish  the  element  coordinate  system.  The  global 
coordinate  system  is  shown  in  Figure  3,  and  for  this  case  local  coordinate  system  for 
each  of  the  elements  coincides  with  the  global  system.  The  ends  of  the  beam  are  simply 
supported  and  restrained  against  rotation  about  the  global  x axis.  Concentrated  loads 
are  applied  at  the  center  of  the  beam,  node  number  4,  in  the  Y and  Z directions.  The 
team  material  properties  and  the  inflation  pressure  are  also  tabulated  in  Figure  3. 

It  is  desirable  to  be  able  to  model  curved  members  with  this  finite  element,  so  to 
establish  its  validity  for  this  class  of  structures,  a circular  arch  was  modeled  as  shown 
in  Figure  4.  This  arch  has  the  form  of  a semicircle  with  a radius  of  120  cm  and  the 
circular  cross  section  has  a radius  of  4 cm.  This  arch  is  modeled  with  eighteen  elements 
of  equal  length.  The  element  members  are  inclused  in  square  symbols  and  the  node  points 
are  denoted  by  the  solid  circular  symbols.  Node  number  twenty  is  a dummy  used  to 
establish  the  local  element  coordinate  system.  The  global  coordinate  system  is  shown 
m Figure  4 and  the  arch  material  properties  and  inflation  pressure  are  tabulated  there 
also  The  ends  of  the  arch  are  simply  supported,  and  displacements  in  the  X direction 
for  all  nodes  are  suppressed.  A concentrated  load  of  unit  magnitude  is  applied  at  node 
number  ten  in  the  negative  Z direction. 

For  both  the  beam  and  the  arch  we  present  the  input  data  and  the  output  results 
from  the  NON F ESA  program  and  a comparison  of  the  finite  element  results  with  the 
results  from  the  exact  solutions  of  the  theory.  The  computer  input  and  results  associated 
with  the  NONFESA  program  are  given  in  Appendices  B and  C for  the  beam  and  arch 
respectively  and  the  comparison  of  the  finite  element  results  with  those  from  the  exact 
solution  are  given  in  Table  1 and  2.  The  input  data  in  Appendices  B and  C are  presented 
as  illustrative  examples  and  no  comments  are  needed.  The  output  from  the  NONFESA 
program  also  appearing  in  Appendices  B and  C can  be  divided  into  two  parts,  the  first 
being  a printout  of  the  input  data  and  some  computed  setup  data  such  as  the 
correspondence  between  degrees  of  freedom  and  equation  numbers.  It  should  be  noted 
that  the  cross-section  radius  and  the  inflation  pressure  are  printed  under  the  Beam 
Geometric  Properties  heading  as  the  Area  X and  Area  Y,  respectively.  Also  note  that 
the  parameter  Inertia  Y under  this  same  heading  is  a negative  number.  As  discussed  above 
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TABLE  2 

Comparison  of  Finite  Element  Resu'ts  with  Those 
from  the  Exact  Solution  for  the  Arch  Check  Case 


Normal  Displacement  Axial  Displacement  Rotation 

cm  cm 


Angular 

Position 

Node 

Finite 

Finite 

Finite 

Degrees 

No. 

Element 

Exact 

Element 

Exact 

Element 

Exact 

-90 

1 

o. 

o. 

o. 

o. 

1.4X10  4 

1 4X10 

-80 

2 

- 3.58X10  ' 

- 3.58X10 

-0.52X10 

-0.53X10 

1.3 

1 2 

-70 

3 

- 5.82 

- 5.81 

-1.57 

-1.59 

0.7 

0.7 

-60 

4 

- 6.39 

- 6.35 

-2.86 

-2.91 

0 06 

0.05 

-50 

5 

- 5.16 

- 5.11 

-4.11 

-4.16 

-0.6 

0.6 

-40 

6 

- 2.19 

- 2.12 

-4.98 

-5.03 

-1.3 

-1.3 

-30 

7 

2.28 

2.39 

-5.19 

-5.24 

-1.8 

-1.8 

-20 

8 

7.89 

8.0 

-4.49 

-4.54 

-2.0 

-2.0 

-10 

9 

13.5 

14.0 

-2.17 

-2.78 

-1.6 

-1.6 

0 

10 

19.2 

19.3 

0. 

0. 

0. 

0. 

Angular 

Position 

Node 

Axial  Force,  N 

Finite 

Moment  N-cm 

Finite 

Stress  N/cm 

Finite 

Degrees 

No.. 

Element 

Exact 

Element 

Exact 

Element 

Exact 

-90 

1 

- 0.525 

- 0.500 

0.0 

00 

-0.020 

-0.019 

-80 

2 

- 0.543 

- 0.545 
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-0  039 

-0.039 

-70 

3 

- 0.572 

- 0.573 

1.31 

1.27 

-0.049 

-0.049 

-60 

4 

- 0.583 

- 0.585 

1.47 

1.43 

-0.052 

-0.052 

-50 

5 

- 0.576 

- 0.578 

1.46 

1.41 

-0.052 

-0  051 

-40 

6 

- 0.552 

- 0.554 

1.29 

1.24 

-0.047 

-0.047 

-30 

7 

- 0.511 

- 0.513 

0.88 

0.836 

-0.037 

-0.037 

-20 

8 

- 0.455 

- 0.456 

0.046 

0.007 

-0.019 

-0.018 

-10 

9 

- 0.385 

- 0.386 

-1.66 

-1.69 
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10 

- 0.347 
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\ 

j 


under  the  description  of  input,,  this  is  used  for  testing  whether  the  beam  element  is  a 
pressure  stabilised  element  for  computation  of  stresses.  The  second  part  of  the  output 
contains  the  computed  results  for  the  beam  stresses,  beam  forces  and  moments,  nodal 
displacements  and  rotations,  and  the  nodal  displacements.  As  is  indicated  on  the  printout, 
for  pressure  stabilized  beam  elements  the  stress  should  he  interpreted  as  a membrane  stress 
with  units  of  force  per  unit  length. 

The  comparison  of  the  finite  element  results  with  those  given  by  the  exact  solution 
of  the  beam  theory  of  reference  1 is  presented  in  Table  1 The  simple  six  element  finite 
element  model  is  in  agreement  with  the  exact  solution  to  all  decimal  places  shown  for 
all  three  parameters,  displacement,  rotation,  and  bending  stress.  The  limit  on  the  number 
of  places  given  in  Trble  1 is  the  output  of  the  NONFESA  program.  The  maximum  possible 
percentage  difference  between  the  finite  element  results  and  tha  exact  solution  is 
then  2 4%.  This  is  good  agreement  and  establishes  the  validity  of  the  element  for  straight 
beams 

The  comparison  of  the  finite  element  results  with  those  given  by  the  exact  solution 
of  the  arch  theory  of  reference  3 is  presented  in  Table  2.  Results  are  presented  for 
half  the  arch  since  the  results  for  the  other  half  can  be  inferred  from  those  presented 
by  the  symmetry  of  the  structure  and  the  loading.  The  finite  element  results  generally 
follow  the  results  from  the  exact  solution  with  differences  of  less  than  5%  and  in  most 
cases  considerably  less  than  5%.  There  are,,  however,  three  cases  where  the  difference 
■s  greater  than  5"o  Two  of  these,  the  axial  displacement  at  node  9 and  the  axial  force 
at  node  10,  do  not  appear  to  have  an  easy  explanation  but  the  third  one,  the  moment 
at  node  8,  appears  to  be  due  to  the  closeness  of  the  node  to  the  point  where  the  moment 
changes  sign  In  this  case,,  the  exact  solution  is  very  close  to  the  value  zero,  and  thus 
any  error  in  the  finite  element  solution  is  magnified.  In  general,  this  comparison  is  good 
and  establishes  the  validity  of  the  finite  element  for  curved  structures. 

CONCLUDING  REMARKS 

The  development  of  a finite  element  for  pressure-stabilized  beams  is  described  along 
with  its  adaptation  to  the  Nonlinear  Finite  Element  Structural  Analysis  program.  This 
provides  a means  for  carrying  out  design  analysis  of  pressure-stabilized  frame-supported 
tents  under  rather  general  loading.  The  validity  of  the  finite  element  and  the  correctness 
of  its  adaptation  to  the  structural  analysis  program  is  established  by  the  comparison  of 
finite  element  results  with  exact  solution  for  a straight  beam  and  a circular  arch., 
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APPENDIX  A 


ELEMENTS  OF  THE  8x8  BENDING  STIFFNESS  MATRIX 


The  elements  of  the  8x8  bending  stiffness  matrix  are  given  here  >n  terms  of  the 
4x4  submatrices  K,  , K,;,  and  K::.  In  presenting  these  elements,  use  is  made  ur 
the  symmetry  of  the  matrices  wherever  oossible 
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NODE 


SYMBOLS 


C 

C,, 

C \ ^ 

D 

E 

G 

h. 

K 

|J 

L 

P 

P 

Q 

S 

S 

|J 

u 

U,,u2 

v,w 

v , t V j ,w , ,w2 

x 


Beam  cross-section  s 
shear  stiffness 
elastic  modulus 
shear  modulus 
bending  stiffness 
axial  stiffness 
torsional  stiffness 

shape  functions  used  in  displacement  approximations 
submatrices  of  the  8x8  bending  stiffness  matrix 
beam  length 
pressure  parameter 
pressure 

vector  of  generalized  forces 
4x4  bending  stiffness  matrix 
elements  of  the  matrix  S 
axial  displacement 
nodal  axial  displacements 
bending  displacements 
nodal  bending  displacements 
axial  coordinate 
vector  defined  by  equation  (9) 
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cross-section  rotations 


nodal  torsional  rotations 
-02,^1  ,^2  nodal  bending  rotations 
< > denotes  differentiation  with  respect  to  x 
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